Ground states of spin-3 Bose-Einstein condensates for conserved magnetization 
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We calculate the ground states and ground state phase diagrams of Bose-Einstein condensates 
of spin-3 atoms under the assumption of conserved magnetization. We especially concentrate on 
the ground states of a 52 Cr condensate. In 52 Cr the magnetic dipole-dipole interaction enables 
magnetization changing collisions, but in a strong magnetic field these are suppressed. In the 
calculation of the phase diagrams we neglect the contribution from the dipole-dipole interaction, 
but discuss its effects at the end of the paper. We show that the ground state of a 52 Cr condensate 
does not seem to depend on whether or not the dipole-dipole interaction is taken into attention. 
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I. INTRODUCTION 



Interest in spin-3 Bose-Einstein condensates (BECs) 
has increased rapidly since the creation of a chromium 
condensate In a chromium condensate the magnetic 
dipole-dipole interaction is much stronger than in alkali- 
metal atom condensates. Therefore two interactions need 
to be considered: the short-range (van der Waals) inter- 
action modeled as a contact interaction, and the long- 
range magnetic dipole-dipole interaction. The latter is 
known to enable conversion of spin angular momentum 
into orbital angular momentum, which makes possible 
the creation of vortex states via spin dynamics [2|,|3j]. The 
spin dynamics driven by the dipole-dipole interaction can 
be prevented by exposing the condensate to a strong ex- 
ternal magnetic field, which suppresses the conversion of 
spin angular momentum into orbital angular momentum 
@- H, IU • The dipole-dipole interaction also affects the ex- 
pansion dynamics of a chromium condensate [5|,[il .allows 
a precise measurement of the scattering length |7|, [8( , and 
opens the way for a new cooling method 9j . One differ- 
ence between 52 Cr and alkali-metal atoms is that in the 
former the nuclear spin is zero. Therefore the quadratic 
Zeeman effect is absent when a chromium condensate is 
placed in a magnetic field. Recently it was shown that 
by using laser fields an effective quadratic Zeeman effect 
can be created also in 52 Cr 10]. 

In experiments performed so far, the condensate has 
been prepared in the S — 3, mj = —3 state, where S 
is the total electronic spin and ms gives the projection 
of S in the direction of the magnetic field. In this state 
the effect of the dipole-dipole interaction is the largest, 
but even then it is weaker than the van der Waals in- 
teraction. Therefore it has been neglected in the studies 
of the ground states of a 52 Cr condensate @, [ll|. This 
simplifies the study of the system, as it allows one to con- 
centrate on the effects of the contact interaction alone. 
Since the contact interaction does not change the value 
of magnetization, the ground states should be calculated 
taking this into attention. For spin-1 condensate this has 
been done in Ref. (l2T |. 

In this paper we calculate the ground states of spin- 
3 condensates under the assumption of conserved mag- 



netization. The ground states of spin-3 condensates 
have been studied earlier, but in the previous studies 
the magnetization has been allowed to change freely 
0, [Tl|, [ll, [T3, We also present the ground-state 

phase diagrams of spin-3 condensates corresponding to 
several fixed values of the magnetization. Finally we dis- 
cuss the special case of 52 Cr and study qualitatively the 
role of the magnetic dipole-dipole interaction. 

This paper is organized as follows. In Sec. |TT]we de- 
rive an expression for the ground-state energy of a spin-3 
condensate. In Sec. Mil we calculate the ground states 
and the ground state phase diagrams. In Sec. IIVI we 
study the effects of the dipole-dipole interaction, analyze 
the ground states of 52 Cr, and discuss the possible ex- 
periments. Section [V] contains a summary of the results 
of the paper. 



II. THE INTERACTION ENERGY 

In this section we derive an expression for the energy 
of a spin-3 condensate. The contact interaction between 
two condensate atoms can be written as U(r — r') = 
S(r — r') X) s =o 9sP s , where V s is the projection operator 

to a two-particle state with total spin s, g s = ^ s 
and a s is the s-wave scattering length in the total spin s 
channel. Because of symmetry, only even terms appear 
in the sum (f6j . We consider a two-particle system and 
denote by Si the spin operator of particle i, that is, Si = 
S (E> ^2, S2 = ii <£> S and S is the spin operator of a spin- 
3 particle. The total spin operator is Stot = Si + S2- 
The possible values for the total spin of two particles are 
now 0, 2, 4, 6 and therefore we get S 2 ot = J2 s =o 2 4 6 s ( s + 
l)V s = 6V 2 + 2OV4 + 42"P 6 {h = I). The sum' of the 
projection operators gives the identity operator, I = 1\ ® 
h = V +V 2 +V 4 +V e . Now Si-S 2 = (S 2 ot -S 2 -S 2 )/2 - 
(S 2 ot - 247)/2 = -12V - 9V 2 - 2P 4 + 9V 6 . By squaring 
this we get (Si ■ S 2 ) 2 = UAV + 81V 2 + 4P 4 + 81V 6 = 
811 + 63"Po — 77P4. The latter expression is obtained by 
using the equation Vq = I — Vo — V 2 — V±. 

The projection operator Vq can be written as Vq = 
|00)(00|, where |00) = \S tot = 0,m Stot = 0). With the 
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help of the expressions for I,Vq, Si • S2, and (Si • S2) 2 
we can write V 2 , V 4 and T 6 as V 2 = - [181 + 2WT + 
7Si • S 2 - (Si ■ S 2 ) 2 ],7> 4 = ^[81/ + 637> - (Si • S 2 ) 2 ] 
and V e = ^[181-30^0 + HSi • S 2 + (Si • S 2 ) 2 ]. Using 
these the interaction term can be written as 



U(r-r') ^6{r-r r )[aI+7f3V +jS r S 2 +S{S 1 -S 2 ) 2 }, (1) 

where a = ~jg 2 +jfg4+i\g6, 7/3 = ffo-f 52+^54-^56, 
7 = 1^(56 - 52) and S = -^ m g 2 - ±g 4 + j^ge- The 
measured values for the scattering lengths of 52 Cr are 
a 2 = —7aB,a 4 = 58ag,aQ = 112(13, while the value of 
do is unknown (as is the Bohr radius) [171] . We write 
a = 47r ^° B a', and define ft, 7', 5' similarly. The scatter- 
ing lengths give a' = 72.19as, 7ft = 42.15as + a , 7' = 
6Ma B ,5' = -0.24a B . 

Next we study the energy of a spin-3 condensate us- 
ing mean-field theory. We write tp = y/n£, where n is 
the particle density and £ is a complex vector with seven 
components and — 1. We call £ a spinor. Previous 
studies have shown that the contribution from the mag- 
netic dipole-dipole interaction is much smaller than that 
of the contact interaction Here we also neglect the 

magnetic dipole-dipole interaction term at first, but dis- 
cuss its effects in Sec. IIV1 If the dipole-dipole interaction 
is neglected the Hamiltonian reads 



which is proportional to 



H = Hn 



dr—ia + CuC 



7(3 (3/x; 3r|00) (00|3//; 3r) 



-7Si w ' ■ S 2rr 



Here H = J dr[h 2 / '(2M)Vip* ■ + Un] with U an ex- 
ternal potential, repeated index is summed, i, j = x, y, z, 
and S iMM ' = S MM / (g) 7 2 , S 2rr / = I x ® S TT /. 

The state |00) appearing in the projection operator P 
can be written in terms of one-particle states as |00) = 
^ [|3 - 3; 33) + |33; 3 - 3) - |3 - 2; 32) - [32; 3 - 2) + 
|3 - 1; 31) + |31; 3 - 1) - |30; 30)]. Using this the energy 
becomes 



E\tb] = 



+ / 3|e| 2 + 7 <S) 2 



(3) 



where (S) = ftS£, O y = (SiSj) = tfSiSjt, and 6 = 
2£ 3 £-3 ~ 2£ 2 £_ 2 + 2£i£_j - C 2 . 

We assume that in the ground state £ is position inde- 
pendent. This is equivalent with the single mode approxi- 
mation, where the orbital state of each spin component is 
taken to be the same. This is a realistic assumption if the 
external potential U is the same for all spin components, 
which is the case in an optical trap. The ground-state 
spinors are then determined by the interaction energy, 



b\ef 



j_ 



dO l 



(4) 



This is obtained by taking the terms in square brackets 
in Eq. ([3]), dropping the constant a and dividing the 
resulting equation by |7|. Following the notation of 
we have defined O 2 = 2 } ■ In the energy < |6| < 1 
and 46 < O 2 < 85.5. The lower bound of O 2 is achieved 
for example when £ = (0,1,0,0,0,0,0) and the upper 
bound when £ = (1, 0, 0, 0, 0, 0, l)/\/2. 

Above we have assumed that there is no external 
magnetic field present. If condensate of spin-3 atoms 
is exposed to a homogeneous magnetic field given by 
B = Be z , to first order in the magnetic field the energy 
changes by Es[ip] = — / drpn(r)(S z ). Here p = g^sB, 
Hb is the Bohr magneton, and g is a constant specific 
to the condensate atoms. We assume that the external 
field is such that it is enough to include this first-order 
contribution only. Since in 52 Cr the nuclear spin is zero, 
for chromium this term gives exactly the energy related 
to the magnetic field. In the presence of an external 
magnetic field the conservation of magnetization has to 
be taken into account. The magnetization is defined as 
M = J drn(r)(S), where in the spin operator H = 1. If 
only the contact interaction is considered, the projection 
of magnetization in the direction of the magnetic field is 
a conserved quantity. Thus now the relevant quantity is 
M z = fdrn(r)(S z ). 

We calculate the ground states under the assumptions 
that the magnetization is fixed and the spinor is position 
independent. The z component of magnetization is then 
determined by m = (S z ) — M z /N, where N is the par- 
ticle number, and Eb[*Ij] = —pNm = —pM z . From this 
on we call m the magnetization. Since (S z ) = m is con- 
stant, we can drop the term proportional to (S z ) 2 from 
the energy, which then becomes 



£ m [V] = b\Q\ 2 + ^r((S x ) 2 + (S y ) 2 ) + dO 2 



(5) 



III. GROUND STATES AND PHASE 
DIAGRAMS 

A. Free magnetization 

We minimize Eq. ([4]) numerically to find out the zero- 
field ground-state spinors. This minimization has been 
performed also in Ref. [ill ] . but using an alternative 
parametrization of the energy. Therefore our figures are 
not directly comparable with those of Ref. [llj . We have, 
however, checked that the zero-field phase diagrams agree 
exactly when the different parametrization is taken into 
account. The paramterization we use makes it easier to 
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see from the phase diagrams how the ground states de- 
pend on b and d. The notation of the ground states 
follows that of Ref. [H. 



|m| 



(a) [7=0, d < 0; 
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FIG. 1: The ground-state phase diagrams if 7 = 0. In (a) for 
d < and in (b) for d > 0. 



Next we briefly discuss the ground-state phases. If 
7 > it is in general favorable to have (S) = 0, whereas 
if 7 < it is energetically advantageous to make | (S) | as 
large as possible. This is why the F and FF phases are 
favored in the figures correponding to 7 < 0. If b, d < 
0, the energy is minimized when |0| and O 2 have their 
maximum values, i.e., |9| = 1, O 2 — 85.5 and the system 
is in A-phase. If d < and b is positive it is favorable 
to have O 2 = 85.5 and |0| = 0. However, these can 
not be achieved simultaneously. If 7 > 0, this leads to 
competition between several different phases, whereas if 
7 < 0, FF is the ground state. If b, d > it is favorable 
to minimize O 2 and |0|. This achieved if the ground 
state is F, for which O 2 = 46 and 9 = 0. If 7 > F is 
replaced by E if d is small, since the latter has (S) = 0. 
Finally, if b < 0, d > 0, it is favorable to set |0| = 1 and 
O 2 =46. These cannot be achieved simultaneously, and 
the ground state is D, for which |9| = 1,0 2 = 48 and 
(S)=0. 



B. Fixed magnetization 

We have minimized the energy shown in Eq. ([5]) us- 
ing m = 0,0.5,1.5 and 2.5 as the values of magneti- 
zation. The ground states are shown in Table |TJ Be- 
cause many of the ground-state spinors resemble the 
zero-field ground states, we have used the same letters 
to label them, but have added the subscript m which 
refers to a fixed magnetization. The ground states de- 
pend on the signs of b, d in a similar fashion than in 
the case of free magnetization. In most phases (S) is 
parallel to the magnetic field, i.e., (S) = (S z )e z . The 
exceptions are F^, FF^, H^, and phases, for which 
(S x ) 2 + (S y ) 2 > 0. These states cannot be ground states 
if the magnetization is free, since then it is favorable to 
rotate the spin so that in the ground state (S) is parallel 
to the magnetic field. An explicit form for F^, FF^ and 
can be obtained with the help of the spin rotation 
matrix i?(ei,e 2) e 3 ) = e ~ MlSs e"' £2S »e" ,E3S «. The value 
of £2 is fixed by the values of d and m, whereas t\ and 63 
are free. For example for the first spinor in Table |T] 
€2 = arccos[— (d+l)m/3d] or £2 = arccos[(d+l)m/3d]+7r. 
These result in different spinors. For FF^ spinor 62 — 
arccos(m/3) and for F^ €2 = arccos(m/2). In Z m and 
Z^ spinors all spin components are populated and the 
exact form of the spinor depends on b and d. 

The energy of a condensate does not depend on the 
sign of m, only on the magnitude. Hence phase dia- 
grams with m = —0.5,-1.5,-2.5 are similar to those 
with m = 0.5, 1.5, 2.5, respectively. The phase diagrams 
for the special case 7 = are shown in Fig.Q] while those 
for positive and negative 7 are shown in Fig. O The (a)- 
(e) diagrams in Fig.[2]correspond to 7 > 0, whereas those 
for 7 < are shown in (f)-(h). If 7 < the diagram with 
free magnetization and that with m = are similar and 
the differences between \m\ = 0.5 and |m| = 1.5 dia- 
grams are very small, so we have omitted the \m\ =0 
and \m\ = 0.5 diagrams. The phase diagrams change at 
m = 2 because D m , E mi and are not possible ground 
states if m > 2, but are replaced by E^. 

In most ground states (S x ) = (S y ) = 0, so the energy 
becomes 



E m [ip] =b\0\ 2 + dO' 2 



(6) 



The ground state of this equation is determined by the 
ratio b/d and the sign of d. If we multiply b and d by a 
positive constant the ground state is unchanged. Equa- 
tion ^ shows that if (d, b) is on the boundary of two 
phases which both have (S x ) = (S y ) = and x is a 
positive number, then also {xd, xb) is on the same bound- 
ary. Therefore the phase boundaries between phases with 
(S x ) = (S y ) = 0, (S z ) = m are straight lines with one 
end at the origin. This reasoning does not work if (S z ) 
depends on b or d, like in B and G phases. 

On a phase boundary the energies of two phases, la- 
beled now by 1 and 2, are equal. Then Eq. Q gives 
b = -d(Of-Of)/(|6i| 2 - |9 2 | 2 ) and by using the values 
for O 2 and |9| 2 listed in Table U equations for some of 
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TABLE I: Here are the ground states. Here m = (S z ) and c (r) is a complex (real) number. R is a spin rotation given 
by -R(ei,£2,£3) = e~" £lSz e~ le2Sy e~ ,,e3Sz . The value of £2 is determined by m. The energy does not change if ms — ► —ms, 
S = —3, —2, . . . , 3, but the sign of magnetization changes. In the second and third column we give the notation used in Refs. 
[IH and 0, respectively. In spinors a,m are such that the terms inside square roots are non-negative and a is a function of 
b, d, and m. 
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the phase boundaries can be obtained. Next we discuss 
the phase diagrams in more detail. 



1. 7 > 

In Fig. [DJa) the boundary between C and H m is at 
d = —0.15, which corresponds to (S z ) 2 = 0.1. In the H m 
spinor m depends on d according to m = — 3d/ (1 + d). 
A and B are separated by b = — 30d, while the B-C 
boundary is b = -43.83d - 98.78d 2 . A-G and H m -G 
boundaries are given by b — 1.5 — 19.5c? and b — —3.20 — 
63.01o! - 76.97o! 2 + 34.82d 4 , respectively. FF and A are 
separated by b — 9 — 4.5d and F and D by b = 4 — 2d. 

In Fig. (Hb) the boundary between fl^ and is lo- 
cated at d « -0.14. 

In Figs.[lJc)-[DJe) F^,^ and F m are divided by a vertical 
line located at d = — \m\/(3 + \m\). This equation is 
valid if |m| > 0.5, for smaller \m\ H m is replaced by Z^. 
The value of a in -f/^ is a — a/(2- d)/5(d+ 1) in the 
first spinor in Table U and a = ^/(l + 4o!)/3(l + d) in 
the second one. These are defined only for those values 
of d for which the term inside the square root is non- 
negative. The boundary between A rn and Z m is given by 
b = — 22. 5d. The Z m -H m boundary depends on m and 
is given approximately by b = (—36.25 + 4.71|m|)d for 
|m| G [1,3]. I m and D rn are separated by 6 = — 30cL 



2. 7 < 



It can be seen that if 7 < the phase diagrams depend 
only weakly on the value of \m\. The ferromagnetic states 
and FF^ are favored since they allow (S) 2 to have 
large values. In Fig. [2ff) F and D are separated by b — 
-2- Ad. FF and F meet at d = 1/7 and FF-A boundary 
is given by b = —9 — 4.5d. In Fig. [^g) D m and 7 m are 
separated by b — — 30d. 



3. 7 = 



The phase diagrams corresponding to 7 = are shown 
in Fig. [1] There are two phase diagrams, in (a) for d < 
and in (b) for d > 0. Since 7 = 0, there are two free 
parameters, the ratio b/\d\ and magnetization m. The 
phase diagram in Fig. [TJa) can be obtained from those 
of Fig. [U by taking the limit d — > —00 with fixed. 
In this limit the term (S^) 2 + (S^) 2 in Eq. ([5]) becomes 
negligible. If 7 ^ 0, the boundary between j4 to and FF m 
is at b = — 4.5<i + 97/I7I. If d — * —00 with b/d fixed, this 
equation gives b/\d\ = 4.5, which is the phase boundary 
in Fig.[TJa). In the same way, by letting d — > 00 with b/d 
fixed, one obtains the phase diagram of Fig. HJb) . 
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FIG. 2: The ground-state phase diagrams. In (a)-(e) we have 7 > and in (f)-(h) 7 < 0. Here b = f3/\~f\ and d = <?/ 1 -y | . 
Notice that in (a)-(e) the scale of d-coordinate changes at d = 0. In (a) the dot-dashed line is determined by the measured 
values of the scattering lengths of 52 Cr and it is located at d = —0.037. It intersects the A-B boundary at b — 1.11 and the 
B-C boundary at b — 1.53, which correspond to do = 9.2a_g and ao = 28.7as, respectively. In all figures except (a) and (f) 
the dashed line means that the phase boundary depends on the value of magnetization, while the solid line denotes that the 
boundary is independent of magnetization. 
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FIG. 3: The ground states of 52 Cr. At the top of the fig- 
ure we give the scale used in Fig. 5 of Ref. [ill ] in order 
to make a comparison with their phase diagram easier. Here 
oldhI^Idh = 1.35 fe— 1.48. The relation between 6 and the un- 
known scattering length ao is given by ao = —42.15 + 46.28 6, 
where ao is given in units of Bohr radius. 



C. Ground states of 52 Cr 



In Fig. [3] we have plotted the ground states of 52 Cr as 
a function of b and \m\. In this figure d = —0.037, which 
is the value relevant for chromium. As explained above, 
phase boundaries between states with (S x ) = (S y ) — 
0, (S z ) = m are of the form b = fed, where k is a constant 
specific to the boundary. If \m\,j > and —0.14 < d < 
one has (S x ) = (S y ) = 0. Therefore Fig. [3] holds true for 
every d' in the interval (—0.14, 0), but with b replaced by 
b' = bd' /dcr, where b is the coordinate used in Fig. [3J 
and dcr = —0.037. If we assume that the uncertainty 
in the values of 02,04, ag is ±10as, then d e [—0.08,0]. 
Thus Fig. [3J is valid for 52 Cr also if a small uncertainty 
in the values of scattering lengths is allowed. 

It is interesting to note that Fig. [3] closely resembles 
the Fig. 5 of Ref. [1JJ. The latter figure gives the 
ground state of 52 Cr as a function of b and external mag- 
netic field, calculated allowing the magnetization to vary 
freely. The similarity is a consequence of the fact that in 
the figure of Ref. [ll| the magnetization of the ground 
state can approximately be written as m = cB, where B 
is the strength of the magnetic field and c is a constant. 
In general, however, the relation between the strength of 
an external magnetic field and the magnetization of the 
corresponding ground state is not of this form. There- 
fore the phase diagrams in Fig. [2] are not identical with 
phase diagrams showing the ground state in some given 
magnetic field as a function of b and d. 



IV. DIPOLE-DIPOLE INTERACTION AND 
THE GROUND STATES OF 52 CR 

A. Dipole-dipole interaction 

Above we have neglected the magnetic dipole-dipole in- 
teraction. If £ is position independent, the dipole-dipole 
interaction energy is 

E d M = — ^— J dvn{v) J dv ^— ^ ' 

(7) 

where (S) = (S)/|(S)|, x = (r - r')/|r - r'|. For 52 Cr 
Cdd = Mc^SeMs) 2 /^, where /j,q is the vacuum permeabil- 
ity, fiB is the Bohr magneton and g e « 2. 

The effect of the dipole-dipole interaction is to change 
the value of 7 and possibly the direction of (S). In most 
experiments the trap is weak in one direction and strong 
in orthogonal directions, which produces a cigar-shaped 
condensate. If (S) is parallel to the long axis of the trap, 
the cloud remains cigar shaped also after the introduc- 
tion of the magnetic dipole-dipole interaction [l^, [20^ . 
Based on Eq. ([7]) one can argue that in a cigar-shaped 
condensate the minimum of Edd is obtained when the 
spin is everywhere parallel to the long axis of the con- 
densate. Thus in the absence of an external magnetic 
field the dipole-dipole interaction aligns the spin with 
the long axis of the condensate [ill ]. 

The assumption of a position independent spinor as 
the ground state is not always realistic. In a 87 Rb con- 
densate of spin-2 atoms the quadratic Zeeman term may 
favor ground states consisting of regions in different spin 
states. The ground state with M z = may, for exam- 
ple, consist of regions where the projection of spin is 
either ms = 2 or mc = —2, and the total magnetiza- 
tion adds up to [18l |, This kind of effect can also be 
caused by the magnetic dipole-dipole interaction. Let us 
assume that the long axis of a cigar-shaped condensate 
and magnetic field are parallel. In most ground states 
(S) is parallel to the magnetic field, which is the direc- 
tion that minimizes the dipole-dipole interaction energy. 
However, for F^, FF^, and states this is not 
true, and therefore the dipole-dipole interaction energy 
is not minimized. The minimum of energy is obtained 
by forming regions where the spin is either along the 
positive or negative z axis. The magnetization of these 
regions has to add up to M z . For example for FF^ state 
the regions would be either in ms = 3 or ms — —3 state. 
Forming these regions costs some kinetic energy, which 
has to be overwhelmed by the energy released from the 
dipole-dipole interaction. 

B. Ground states of 52 Cr 

The scattering lengths of 52 Cr indicate that in the ab- 
sence of dipole-dipole interaction and for free magnetiza- 
tion the ground state is A, B, or C . If the magnetization 
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is fixed and \m\ is larger than 0.5, the ground state is ei- 
ther A m , G m , H mi or Z m . The difference between these 
states is easy to see experimentally, since they have dif- 
ferent spin components populated. 

The introduction of the dipole-dipole interaction may 
change the ground states. For chromium the value of 
7 is positive and 5 is negative. We assume that the 
change in 7 following from the dipole-dipole interaction 
is such that 7 remains positive. A change in 7 shifts the 
point (-A- , ) characterizing chromium radially toward 
or away from the origin. It, however, cannot cross the 
origin. 

If magnetization is free, one can see from Fig.[5]that the 
H m phase is reached if S/-f < —0.15 and b is large enough. 
For 52 Cr these conditions are fulfilled if the change A7 = 
— IA7I induced by the dipole-dipole interaction is such 
that I A7I/7 > 0.76 and b > 1.18. The latter is equivalent 
with ao > 12.Qas- 

If magnetization is fixed, the ground state is either 
unchanged or becomes H**, FF^, or as the value of 
7 changes. The amount of change in 7 required to reach 
H m depends on the value of magnetization. If |m| < 
0.5, the boundary between and is at d m —0.14, 
which requires | A7I/7 > 0.74. For \m\ > 0.5 the H m - 
boundary is at d = — \m\/(3 + \m\), which shows 
that IA7I/7 > 0.96 — 0.11/|m| > 0.74 is required in order 
to reach the state. Also the minimum value of b for 
this to be possible is a function of |m|, but if b > 1.18, or 
ao > 12.6as, H m can be reached, regardless of the value 
of m. 

The change in 7 caused by the magnetic dipole-dipole 
interaction can be estimated by using the analytical re- 
sults for the dipole-dipole and contact interaction en- 
ergies obtained in [l|| [2(J. These results hold in the 
Thomas-Fermi limit, where the kinetic energy is negli- 
gibly small compared with the interaction and trapping 
energies. We assume that (S) = me z , the trapping po- 
tential is cylindrically symmetric, U(r) = M[uj 2 (x 2 + 
y 2 ) + uj 2 z 2 ]/2, and that lo z /lo x < 1, so that the trap is 
cigarlike. In Fig. [5] we plot the ratio A7/7 as a function 
of the trap aspect ratio u) z /u> x . The minimum value is 
A7/7 w —0.26. A similar result was found in [ll[ using 
a Gaussian ansatz for the particle density. In addition to 
the dipole-dipole interaction, an uncertainty in the val- 
ues of the scattering lengths may produce a negative A7. 
If the values of 02, 04, and ciq are known within ±10as, 
the maximum change is IA7I/7 = 0.17. Thus the change 
in 7 resulting from the dipole-dipole interaction and a 
possible uncertainty in the scattering lengths does not 
appear to be big enough to change the ground state of a 
chromium condensate. 



C. Experiments with Cr 

In a high magnetic field the coherent spin relaxation 
via magnetic dipole-dipole interaction is suppressed. 
This occurs when the energy spacing of the neighboring 
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FIG. 4: The change in 7 induced by the dipole-dipole inter- 
action. Solid line is for |m| = 3 and dashed line for m = 0. 
These give the maximum and minimum for IA7I, respectively. 
If m = 0, the ground state is independent of A7. 



Zeeman levels, g e HBB, is much larger than the chemical 
potential fi. For a chromium condensate of 10 5 particles 
these energies are equal when B ~ 2 mG. Although co- 
herent spin relaxation is prevented, spin relaxation can 
still occur via incoherent processes. In these the energy 
released in spin relaxation is so large that the colliding 
particles escape from the trap. For 52 Cr this process 
has been studied both theoretically and experimentally 
in Ref. [2ll| . For example for 777,5 = 3 state at a magnetic 
field of 1 G the rate constant pdr characterizing inelas- 
tic dipolar relaxation was found to be (3dr — 4 x 10~ 12 
cm 3 s _1 [2l|. Thus for a typical central density 77 = 10 14 
cm -3 the 777s = 3 state decays in a few ms. This is a 
short time compared with the typical timescales of 52 Cr 
experiments, which are on the order of tens of ms @, H2|. 
One way to suppress the incoherent dipolar relaxation is 
to reduce the condensate density. This can be done by 
loading the condensate to an optical lattice of shallow 
traps [23[. If one reduces the density, e.g., by two orders 
of magnitude, the new decay time is on the order of hun- 
dreds of milliseconds, which should be long enough for 
experiments. 

From experimental point of view it is important to 
be able to produce different spinors as initial states. 
This can be done using rapid adiabatic passage and con- 
trolled Landau-Zener crossing techniques, as described 
in Ref. [24j. For spin-1 and spin- 2 87 Rb condensates 
this has been done in Ref. |25|. For | ttt. | > 0.5, 
the ground sta te of 52 Cr is very likely of the form 
(0, v/(3 -(- 777)/5, 0, 0, 0, 0, y/(2-m)/5). If a condensate 
is initially prepared in this state, it should not show 
any coherent spin dynamics. Because of incoherent pro- 
cesses the populaton of the 7775 = 2 component decreases, 
while that of 7775 = —3 component is almost unchanged. 
This leads to a decrease of magnetization. If initially 
777 < —0.5, at every later moment the condensate is in 
the instantaneous ground state of spin with magnetiza- 
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tion m < —0.5. 



V. CONCLUSIONS 

We have calculated the ground states and the ground- 
state phase diagrams of spin-3 BECs using mean-field 
theory and single-mode approximation. We have as- 
sumed that the projection of the magnetization in the 
direction of the external magnetic field is a conserved 
quantity. We have presented the ground-state phase dia- 
grams pertaining to several values of magnetization. The 
phase diagrams are classified by the value of 7, which is 
determined by the scattering lengths of different scatter- 
ing channels. We found out that the phase diagrams with 
7 < depend only weakly on the value of magnetization, 
whereas if 7 > the dependence is stronger. The phase 
diagrams for a 7 = condensate are obtained in a certain 
limit from the 7 > or the 7 < diagrams. 

We have plotted the ground states of 52 Cr as a function 
of magnetization and the unknown scattering length a®. 
We found out that this resembles the phase diagram for 
a chromium condensate in a weak magnetic field, where 
the magnetization can vary freely. This follows from the 
fact that in the latter the magnetization of the ground 
state depends approximately linearly on the strength of 



the magnetic field. 

We have also studied how the magnetic dipole-dipole 
interaction affects the ground states. In the absence of 
the dipole-dipole interaction it is possible to have ground 
states where the magnetization is not parallel to the 
magnetic field. We have showed that, instead of these 
states, the dipole-dipole interaction may favor ground 
states where the condensate has broken up into regions 
having the magnetization parallel to the magnetic field. 
We have illustrated that the ground state of 52 Cr does 
not appear to depend on whether or not the contribution 
from the dipole-dipole interaction is included. Finally we 
have discussed the experimental realization of spinor 52 Cr 
condensates, pointing out that due to fast incoherent spin 
relaxation the condensate density should be decreased in 
order to make the lifetime of condensate long enough for 
experiments. 
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